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The key feature of the Fe-based superconductors is their quasi 2D multiband Fermi surface. 
By relating the problem to a negative U Hubbard model and its superconducting ground state, 
we show that the defining instability of such a Fermi surface is the valley density- wave (VDW), 
a combined spin/charge density- wave at the wavevector connecting the electron and hole valleys. 
As the valley parameters change by doping or pressure, the fictitious superconductor experiences 
"Zeeman splitting" , eventually going into a non-uniform "Fulde-Ferrell-Larkin-Ovchinikov" (FFLO) 
state, an itinerant and often incommensurate VDW of the real world, characterized by the metallic 
conductivity from the ungapped remnants of the Fermi surface. When "Zeeman splitting" exceeds 
the "Chandrasekhar-Clogston" limit, the "FFLO" state disappears, and the VDW is destabilized. 
Near this point, the VDW fluctuations and interband pair repulsion are essential ingredients of 
high- Tc superconductivity in Fe-pnictides. 



INTRODUCTION 

Recently, the superconductivity below 7 K in LaOFeP 
[H led to the discovery of high 3^ ^ 26 K in its doped sib- 
ling LaOi_3:Fa;FeAs {x > 0.1) (4I. Even higher T^'s were 
found by replacing La with other rare-earths (RE) , up to 
the current record of Tc = 55 K 0] . These are the first 
non-cuprate superconductors exhibiting such high Tc's 
and their discovery has touched off a storm of activity 

a. 

In this paper, we introduce a new element into the 
theoretical debate by considering a unified model of spin 
density-wave, orbital density-wave, structural deforma- 
tion and superconductivity in Fe-pnictides. The model 
is simple but it contains the necessary physical features. 
The essential ingredient are electron and hole pockets 
(valleys) of the quasi two-dimensional (2D) multiply- 
connected Fermi surface (FS) [1, @, 0]. To extract the 
basic physics we consider spinless electrons first, and 
only a single electron and a single hole band with iden- 
tical band parameters. We then show that this model 
can be related to a 2D negative U Hubbard model, 
the ground state of which is known exactly - it is a 
superconductor [sj. In real FeAs materials, this ficti- 
tious superconductivity translates into a fully gapped 
valley density- wave (VDW), a unified state representing 
a combination of spin, charge and orbital density-waves 
(SDW/CDW/ODW) at the commensurate wavevector 
M connecting the two valleys. Next, we introduce two 
different "chemical potentials", /x*^ 7^ /x'' for the elec- 
tron and the hole valleys - this describes the effect of 
doping the parent iron-pnictide compounds and corre- 
sponds to the external Zeeman splitting in our ficti- 
tious negative U Hubbard model. As 6fj, — fi"^ — /x'' in- 
creases, so does this Zeeman splitting, and eventually our 
fictitious superconducting state approaches to and ex- 
ceeds the "Chandrasekhar-Clogston" limit, giving way to 
a non-uniform Fulde-Ferrell-Larkin-Ovchinikov (FFLO) 
ground state at an incommensurate wavevector q, where 



|q| is set by 6kp = kp — kp, and thus by doping x. This 
"FFLO state" is nothing but an incommensurate (IC) 
VDW at the wavevector M -I- q. Finally, as Skp (x) ex- 
ceeds certain critical value Skc (xc), the "superconduct- 
ing" state is completely destroyed and so is the VDW 
in a true material. However, for Skp above but near 
Skc, we consider strong "superconducting" fluctuations 
and find that these VDW fiuctuations can induce real 
superconductivity in Fe-pnictides (see Fig. [1]). In princi- 
ple, one could avoid the mapping to the negative-U Hub- 
bard model and argue that the VDW instability in pnic- 
tides occurs for the same reasons as the SDW instability 
found in, say, Cr [9]. We find, however, that our 'ficti- 
tious superconductivity' description is more appropriate 
to pnictides not only due to its illustrative purposes, but 
also because it allows us to extend the analogy to the 
"FFLO" state and multiband "SC", i.e., VDW. Furthe- 
more, it provides a natural venue, by using the known 
near -rig ourous results on two-dimensional superconduc- 
tors [10], to highlight the crucial role played in real super- 
conductivity by the interband pair-scattering processes, 
as dicussed at length below. 

The appearance of an SDW in pnictides along with a 
structural transition, which we argue to be a signature 
CDW coupled to phonons, has been established early on 
in the so-called 1111 family 



111, and the universal pres- 



ence of these orders in other families of pnictides and 
related materials (122, 11, etc.) has been confirmed by 
many authors 13, l3| ■ Both the magnetic and the struc- 



tural order set in at nearly identical temperatures, the 
experimental fact that has motivated us to model this 
problem as a one major, "mother" instability (VDW) 
driven by a large energy scale, which is then split into sev- 
eral stages (say, CDW, followed by SDW, and eventually 
ODW order) by interaction terms considerably smaller in 
magnitude. A main feature of our fictitious FFLO state 
are the ungapped portions of the reconstructed Fermi 
surface(s) (FS's). These FS's have been observed in pnic- 
tides directly [ij] as well as indirectly through their sig- 
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natures in metallic resistivity ^iB*! and recently detected 
incommensurability of the SDW order [l^. 



AN IDEALIZED MODEL OF A VALLEY 
DENSITY- WAVE AND AN FFLO STATE IN A 
FICTITIOUS SUPERCONDUCTOR 

The band structure of Fe-pnictides 0, 0, 0] can be 
parametrized by the five-orbital tight-binding model 
[13, [3] ■ The key feature is the multiband nature of the 
FS, possessing both hole and electron sections. We work 
with the properly defined Fe-pnictide unit cell which con- 
tains two of Fe and two pnictide (As or such) atoms per 
unit cell. The basic physics is captured by the Hamilto- 
nian describing two hole (c") and two electron [d^] bands 
centered at the F and M points of the 2D Brillouin zone 
(BZ), respectively, 

H = Hq + Hint , (1) 
rr _ ,(")^(")t J") , ,(/3) j(/3)t j(/3) /9^ 



k,(T,/3 



H,nt - ^ J d^rd^r'V{r,r')n{r)n{r') 



(3) 



where ct, cr' ==T)i and a, (3 are the spin and band la- 
bels, respectively, e^"''^'' is the hole (electron) disper- 
sion near the FS, V{r,r') is the effective interaction and 
n{r) = E>J(r)^.(r), with ^,(r) = Ek,o^>L"Hr) + 

Ek,/3rfS</'lf^(i-)- and (l>i^\r) are the Bloch 

wavefunctions of hole (electron) bands. 

For simplicity, ^ includes only the screened density- 
density repulsion; its form becomes more complex if we 
integrate out the bands away from the Fermi level Ep, 
generating additional interactions in the spin and inter- 
band (orbital) channels. Furthermore, we could equally 
well start from the minimal tight-binding representation 
of Ref. [3] and introduce the interaction term in the 
Wannier representation as 



Hi7it ^^Ud^nli- Jb 



und E] ^di 



+ (•••), (4) 



where ndi and Sdi are the total particle number and spin 
in Wannier d-orbitals of iron. Ud describes the overall 
Hubbard-like repulsion on iron sites while Jnund signi- 
fies the intra d-orbital Hund coupling. In addition, there 
are numerous intra d-orbital interactions, as well as var- 
ious similar terms for p-orbitals on pnictide sites, all 
contributing to (• • • ) [l9|. However, all such interaction 
terms feed into the generic classes of quartic vertices near 
the FS which are generated by Hint ©; only the precise 
numerical values of various vertices are affected. In par- 
ticular, as long as the influence of Jnund (and (• • • ) terms) 
is relatively small and one is in the weak-to-intermediate 
coupling regime argued to be relavant to pnictides [l8| . 



the overall numerical hierarchy of energy scales defined 
by these various classes of vertices remains intact, as de- 
scussed below. Finally, we further simplify the problem 
by exploiting the fact that all electron (hole) bands have 
Ep near their bottom (top) and their Fermi wavevec- 
tors kp^s are ^ M. This allows us to restrict our at- 
tention to the first BZ and take continuum limit, with 
V{r, r') ^ V(r-r'). 

Hint Q generates three classes of vertices: i) the in- 
traband {c'cc'c and d^dd'^d), ii) the interband [c'cd^d), 
and iii) the mixed {d'' cc^d + h.c. and dc'^ d + h.c.). All 
arise from Hint 5 J2q ^q'^q"--q, where 



EC: 



{aa'} (Q)t (a') 
k+q,k^k+q,(T^k,o- 



EC, 

kcr/3/3 



(1313'} ,(/3)t ,(/3') 
k+q,k"k+q,(T"k,(T 



/ . 'k+q,k^k-|-q,o-"k,(T ^ n.c. 

k(TQ/3 



(5) 



is the Fourier transform (FT) of y(r — r'), and 
CSr/' = / d^re'('''-'^)>L"^*(r)v.L?''(r), 



''k,k' 



d^re^C^'-'')--^^ *(r)0(r)(r), (6) 



<t = J d^re^('''-''^-VL"^*(r)<^L'^^(r). 

The following should be kept in mind about these three 
classes of vortices: first, all exhibit considerable varia- 
tion as one moves around the FS. This is the consequence 
of significant variations in the orbital content of various 
bands in different portion of the BZ. Second, we find that, 
generically, the intra and the interband vertices are com- 
parable in magnitude while the mixed ones are notably 
smaller. This remains true regardless of whether we use 
the interaction ([3]), jl]) or some other related form as long 
as Jnund is not dominating the physics and (• • • ) dH) are 
relatively small. 

To illustrate the latter claim, we include the Hund's 
coupling to the interaction terms and compare the ver- 
tices with and without it. For example, the hole intra- 
band vertex 



' k , k' , q ''k+ q , (T "^k' - q , CT ' '^k' cr ' '^ko- 

acquires strength 

'^k.k'.q ^ I'^q 1^ 4JHund; '5k-|-q,k'>k'-q,k' ^ 

1 J ^(aa) ^{aa) , r.^ 

2>^Hundl,k+q,k''5k'-q,k' \° ) 

and therefore the Hund's coupling effectively only adds 
up to the Coulomb potential in this scattering channel. 
The same is true for other intraband scattering processes, 
the second mixed term G2 , and most importantly for the 
interband scattering vertices 



(7) 



<k?q - (^q + ^Hund) cr+;!kC, 



[aa] Jf3l3) 



k+q.k^'k'-q.k' 



2 '^HundVk+q^k'/k.k'-q' 



(9) 
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where the second term is neghgible for small q. The only 
interaction vertex that is more affected by Jnund than 
the others is the first mixed term 

•^l/k.k'.q - l'^q+ 4'^Hundj7k+q,k7k',k'-q + 

2>^Hundi,k+q,k'^k'-q,k" V-^"^ 

due to the relative size of C's and 7's. However, as long 
as the screened Coulomb potential is the strongest inter- 
action (i.e., Jnund ^ here) the changes to vertices due 
to the other sources of scattering (such as those in Eq. 
(jl])) will be only quantitative in nature. 

We now observe that the shapes of different sections 
of the FS (Fig. [T]) resemble each other to a reasonable 
degree. Furthermore, various masses are also roughly 
similar [13, [3 ■ Thus, to make theoretical progress, it is 
useful to first assume that all electron and hole bands are 
equal - - e^^,^ - ^^^m - <^k- ^t" 

ter making the particle-hole (p-h) transformation — » 
^iTct' \ ~^ ''^''■k'-^T' Hamiltonian ([T|) becomes: 

^su(8)-E^k4'K'*K+^^'-' (11) 

k,cr/i 

where *(^)t = {e\,e\,h\,h\). I gnoring the effects of 
bands away from the FS, the kinetic part of i?su(8) pT|) 
has an exact SU(8) symmetry involving orbital (both elec- 
tron and hole) and spin degrees of freedom, /i and cr, re- 
spectively. This symmetry can be used to classify various 
vertices in - ultimately generated by Hint © which 
itself has only U(l)chargcXSU(2)spin symmetry - and an- 
alyze various symmetry-breaking patternSj_starting with 
SUi^%) SU{A) X S'[/(4), as discussed in 

To illustrate the physics, we focus first on the minimal 
model: SU(8) SU(2). This leaves one with only two 
fermion flavors, e and h. Note that spin SU(2) symmetry 
is suppressed, but the orbital electron-hole symmetry re- 
mains, as it is essential for this problem. We now obtain: 

^^su(2) = ^k[e],ek -I- ^'^k] + W d^rd^r' x 

k 

[U'ir - r')n^{r)n'{r') + t/''(r - r')n''(r)n''(r') 
-2W{r ^ r')n^(r)n^r')] + H^^^^^f , (12) 

where U'^^^{y) and W{r) are the Fourier transforms of 

V^.-^.'&&hs and F^.k^iC^Lc^^Hk', k))ps, re- 
spectively. (•••)fs indicates the k- and k'-dependence 
was replaced by the average over the FS - this is justi- 
fied later. The SU(2) symmetry imphes C/^(r) = U^{y). 
Finally, the general H™^^°'^ contains smaller mixed ver- 
tices and eventually plays a prominent role in our theory; 
however, we initially — but only temporarily — set it to 
zero, H^^^f ^ 0. 

The intraband scattering VF(r) has a minus sign in 
front of it. This is the result of the p-h transformation 




FIG. 1: (Color online) Phase diagram of Fe-pnictides, depict- 
ing the evolution of our fictitious superconductor from the 
fully gapped VDW insulator to the "FFLO superconductor" - 
a partially gapped metallic VDW - to the real superconductor 
(SC) under the influence of the "Zeeman splitting" Sn (doping 
or pressure). Red dot symbolizes the parent compounds and 
the regime below it might be physically inaccessible. Insets: 
FS of a) the normal state in the folded {F ^ M) BZ b) 
the VDW metal (computed with the interband interaction set 
to unity) - this is the C4 version of c) the continuum FFLO 
state [2J]. The remaining states are fully gapped. 



and indicates that, us having started with a (screened) 
Coulomb repulsion between the original electrons, the e 
and h flavors now mutually attract. Consequently, at low 
energies, -ffsu(2) (without H^^°^) is equivalent to the 
negative U Hubbard model at low (or high) x, assum- 
ing that W{r) is short-ranged and ek can be represented 
in the effective mass approximation {W{r) W6{r) 
and Ck k^/2m). Both assumptions should be valid 
since kp <C M. The ground state is a "superconduc- 
tor", with an anomalous correlator (e^h^) ^ 0, where 
e^{h^) creates a "spin up (down)" fermion (/|). For 
T > Tc ^ EFexp{-l/N{0)W), the system is in its nor- 
mal state (see phase diagram in Fig. [1]). At T < Tc one 
enters a broken symmetry state, with "off-diagonal" long 
range order and gapped fermions. 

Of course, the above "superconductivity" is useful but 
entirely fictitious mathematical construct, resulting from 
the p-h transformation used to enhance the symmetry 
of our basic model for iron pnictides. Still, what is 
this "superconductivity" in the real world? By retrac- 
ing our steps and undoing the p-h transformation, the 
"off-diagonal" order in (/kf/lk^) (sjk'^-k) translates 
into the diagonal density-wave, ('^k_|_]v['^k) 7^ 0, connect- 
ing electrons from two pockets of the FS separated by 
M — (tt, 7r); a valley density-wave (VDW). Note that 
the VDW describes both spin and charge/orbital density- 
waves (SDW/CDW/ODW). With the spin SU(2) sym- 
metry suppressed in our minimal model, one cannot - 
and should not - distinguish between the two. We iden- 
tify the above VDW (SDW/CDW/ODW combination) 
formation as the physical mechanism driving the density- 
wave orderings observed in numerous experiments. 
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We can pursue this VDW- "superconductor" analogy 
a bit further: in real FeAs materials, the electron-hole 
pockets are not identical, the main difference being their 
distinct kp^s. In our fictitious superconductor, this trans- 
lates to different "chemical potentials", fi'^ ^ fi^ for 
the electron and the hole valleys. This is nothing but 
the external Zeeman splitting in a fictitious negative U 
Hubbard model. As 6fi = fi'^ — ^jP- grows, the "su- 
perconducting" state approaches the "Chandrasekhar- 
Clogston" limit, giving way to a non-uniform Fulde- 
Ferrell-Larkin-Ovchinikov (FFLO) ground state at an in- 
commensurate wavevector q, where |q| is set by 5kp = 
kp ~ kp. This "FFLO state" is just an incommensurate 
(IC) VDW at the wavevector M -I- q. Finally, as 5kp (x) 
exceeds certain critical value 6kc {xc), the "superconduct- 
ing" state is destroyed and so is the VDW (SDW/CDW) 
in a real FeAs system (Fig. [T]) ^22.] . 



VALLEY DENSITY- WAVE AND INTERBAND 
SUPERCONDUCTIVITY 

The above "superconductor" analysis dealt with an 
idealized model but its main conclusions apply to the real 
Fe-pnictides: i) the dominant instability is the VDW at 
wavevector M, a unified spatially- modulated state man- 
ifested through the combined SDW/CDW/ODW and a 
structural transformation [lljl, the details of which de- 
pend on non-universal features of individual materials; ii) 
since hole and electron valleys are not identical, the VDW 
is the p-h analog of a fictitious FFLO state, resulting al- 
most always in portions of the FS which are not gapped 
(Fig. [I]). Consequently, the SDW/CDW/ODW's in FeAs 
are highly itinerant and coexist with finite density of nor- 
mal charge carriers, exhibiting metallic conductivity [l5| ; 
and, finally, iii) the Hamiltonian -ffsu(2) P2p without the 
mixed vertices (i.e., with H^^^'^ 0) contains only three 
basic ground states: fictitious uniform and non-uniform 
FFLO "superconductor" (C and IC VDW) and the nor- 
mal state (Fig. [T]). Thus, if purely electronic interactions 
are to have a prominent role in generating Fe-based su- 
perconductivity of the real world, this effect must arise 
from H^^^f- This is an important result, and an added 
benefit of our transforming the original problem into a 
fictitious superconductor [10|. 

With this last point in mind, we now restore these 
mixed vertices to investigate the real superconductivity 
in our "superconductor" model. It is beneficial at this 
stage to add extra two fiavors to the elementary SU (2) 
model and demand an additional global SU{2) isospin 
symmetry with respect to these new flavors - this isospin 
degree of freedom is completely inert and can be thought 
of as either the real spin or an additional orbital index. 
Its role is purely mathematical as it couches the following 
analysis in the language most easily translated to the ulti- 
mate realistic description of pnictides (23| . Furthermore, 



viewing this isospin as simply the real spin is useful since, 
assuming a reasonable degree of total spin conservation 
in pnictides, the additional SU(2) global symmetry limits 
the number of terms in H^"^''^ that need to be consid- 
ered. The mixed vertices allowed by this extension of our 
model are mixed scattering Gi , and Josephson-type term 
G2 which, in absence of nonlocal interactions, has to be 
in a spin-singlet channel 



mixed 



Gic\d„dl,c„ + 

\G2{cTcUla){^'d-a'd,.) + h.c. . (13) 



The relation of these vertices to the screened Coulomb in- 
teraction is given in Eq. (|10p . We assume that the corre- 
sponding coupling constants are in the regime G\ , G2 
W,Uh,Ue, as will be justified momentarily. 

These preliminaries in place, we are ready to answer 
the key question: what is the effect of finite Gi, G2 on the 
previous analysis? Analyzing corrections in the pertur- 
bation theory to our four point vertices, we find contribu- 
tions to the processes which are dependent on whether 
the incoming, and outgoing, spins are parallel or not. 
This occurs due to the fact that incoming or outgoing 
states of mixed term G2 as well as the intraband interac- 
tion U , are spin singlets of either holes or electrons, Eq. 
P3p . The interband scattering term is therefore conve- 
niently split into two pieces 



Wdlch 



W'dl 



.ciCr. 



d^, + W"dlcUada, (14) 



with the bare values for both coupling being identical to 
W . The first mixed term Gi is split in identical fashion, 
while the intraband scattering and the G2 mixed term are 
required to scatter particles with incoming opposite spins 
and therefore do not need to undergo the same separation 
(or equivalently, G2 = 0, etc.). At the lowest order, we 
find that different types of vertices receive the following 
corrections in the perturbation theory: 

,A, ,A, 

■■uj'pp -oj' 
„ wA. ^ , . .A. 
92 - 2.g2.9c/ ln(- 



5(7 (^) 
52 (t^) 



2g2g'wH-)lh 



,A 



2g2g'^ln{-)l^-2g2g'{\n 



9wi^) 
9wi^) 
5i(^) 
9i{^) 



9w 



9w 



i9'wrH-] 

-) 

A 



-) 

A. 



ph 



.92ln(- 



UJ 



5i'"(5i)'ln( 



A. 



ph' 



A. 



(.9'i')^ln(-) 



ph 



.A, 



9lHt±)^^ + 29W^H-)lh' 



(15) 



where gu, 52, 9'wi-3w^ 9i and g'{ arc just the vertices 
U{= Uh = Ue) @, G2, W, W", G; and G'/, re- 
spectively, measured in units of inverse density-of-states 
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(DOS) at the Fermi level. The logarithmic divergences in 
(fT5)l arise from two sources: first, the standard Cooper 
pairing instability in the particle-particle (pp) channel 
and, second, the perfect nesting of the hole and electron 
bands in the particle- hole (ph) channel, i.e., our ficti- 
tious "Cooper" instability. Finally, In(^)^^''^ denotes a 
crossing (vertex) diagram in the p-h channel — it strictly 
diverges only in the kp'^ /M limit and is otherwise 
finite. Needless to say, there are many additional terms 
that contribute to various vertices at the leading order in 
perturbation theory. However, all such terms are finite 
in the low energy limit and are omitted from (jlSp . 



The third and fourth lines of (fT5|) are just the mathe- 
matical shorthand for our earlier discussion: under renor- 
malization, the coupling constants g\y and g'^ keep grow- 
ing, ultimately generating the VDW instability, driven 
by W (the short-ranged attraction of our fictitious Hub- 
bard model). We notice, however, that G2 enhances the 
growth of W' , thereby giving slight edge to SDW (spin- 
triplet) over a CDW (spin-singlet). The first and last two 
lines tell us that U and Gi's do not interfere: the intra- 
band repulsion gu^ initially large, is rapidly renormalized 
downwards, toward the Fermi liquid behavior, gis are 
typically small to begin with and are also driven down; in 
practice, they can be set to zero. The interesting physics 
is reserved for G2. The growth of gvi^'s fuels the growth 
of 52 and thus the mixed vertex describing the resonant 
pair scattering between the hole and electron bands — 
i.e., the "Josephson" interband vertex c^dc'^d + h.c. — 
becomes strongly enhanced as one approaches the VDW 
(SDW/CDW/ODW) instability. However, since typi- 
cally G2 W (see below), gw wins, resulting in the 
VDW order. Once the VDW is formed, the fermions are 
gapped and the singular behavior disappears, and with 
it any additional enhancement of 6*2. 



The situation changes, however, when the differences 
in size between h and e bands are included, i.e., when the 
"Zeeman splitting" is turned on, by doping or pressure in 
real FeAs (Fig. [1]) . This cuts off the fictitious "Cooper" 
divergence, resulting in our "FFLO" state and eventual 
disappearance of VDW. In this case, the portions or even 
all of the FS is still available for the true Cooper pairing 
and the real superconductivity becomes a viable option. 
The remarkable feature of interband pair resonance is 
that it can produce real superconductivity irrespective 
of its sign [2J|. Thus, strongly enhanced G2 can take 
advantage of the real-world Cooper singularity - which 
is always present - and amplify a preexisting intraband 
superconducting instability or generate one entirely on 
its own. We will revisit this point shortly. 



VALLEY DENSITY- WAVE AND 
SUPERCONDUCTIVITY IN REAL IRON 
PNICTIDES 

This is as far as we can go within the idealized pic- 
ture: we now must face up to the complexities of the 
real materials. First, there are four, two h and two e, 
bands which deviate from an ideal parabolic shape and 
whose kp's are different and, second, all vertices - in- 
tra, interband and mixed - have considerable structure 
as one moves over different portions of the FS. The latter 
is an important point and reflects a fundamental feature 
of FeAs: all five d-orbitals need to be included in realistic 
calculations and various two- and three-orbitals models 
will fall short in addressing the phenomenology of real 
materials. We find that, no more than a single orbital, 
c^222_2;2_y2, can be dropped without a major qualitative 
disruption of the character of the electronic states at the 
FS; thus a four-orbital description is the absolute mini- 
mum. Finally, the lattice effects produce modifications to 
our continuum picture which need to be addressed 22 1. 

We use the full 8-1-8-band tight-binding model to 

find the electron (y'^"'') ^^'^ iiole {4>'^^) wavef unctions. 
This model yields C's and 7's (O shown in Fig. [2] For 
a fixed k, a given ( varies as k' goes around the FS. At 
k = k', the normalization of wave- functions sets ^ to 
1. As k and k' move apart, so does ( decreases until it 
reaches its minimum at k' = — k. Based on the symmetry 
properties of the atomic orbitals, one obtains {hp <^ M) 



'5k,-k 



where is the amplitude of atomic orbital /i in a 
hole state (a) or, equivalently, an electron state (/3). 
Each orbital's contribution is determined by its in-plane 
parity (i.e., sign change under {x,y,z) — > {—x,—y,z)); 
even/odd orbitals contribute with -I-/—, or vice versa. 
Our model uses orbitals of different parity, and conse- 
quently, is bound between —1 and 1, the precise 
value depending on the amount of mixing of even and 
odd orbitals within a state. For example, compare C^'^) 
and C^''^^ (Fig. [21). Since both hole bands have a signif- 
icant contribution of d^z/yz atomic orbitals (odd), both 
are similarly shaped. The 4-fold repetitive structure is 
due to the C4 lattice symmetry. However, the minimum 
values (k = — k') are different: CC*^) nearly reaches — 1, 
whereas ('C'^) > —0.6. This reflects the fact that the 
outer hole band possess a significant overlap with dxy 
(even) orbital state, while the inner hole band is almost 
entirely made of odd bands. In a more limited model, 
where only bands of a certain parity are kept |25| , a topo- 
logical 'Berry phase winding' can be defined for each sec- 
tion of the FS [2^. Depending on this 'winding', Ck,-k 
would have to be either -1-1 or —1, and the consequences 
of the latter would include a suppression of an s-wave 



)|2^ 



(16) 



fi^even 



fi^odd 
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VDW in a favor of a p-wave one. Within our model, this 
notion of topology is absent. 



a) 




FIG. 2: (Color online) a) inter/intraband Ck^k^"' of Eq. © 
(in fcF ^ M limit) for the inner and outer hole bands, and 
the electron bands as k and k' go around respective FSs. b) 
mixed 7^'^^''' © s-re clearly smaller than ("'s. 



Next, the above form factors C's and 7's are used to 
compute all the interaction vertices (intra, interband and 
mixed) stemming from ([3]) along different sections of the 
FS (/il, /i2, el, and e2) and to extract the corresponding 
coupling constants in the C4 "angular momentum" chan- 
nels, s, p, and d. The results, normalized by the overall 
strength of the screened Coulomb interaction in ([3|), are 







U 






W 






(hi) 


(h2) (ol) 


(o2) 


(hl.ol) 


(hl,o2) 


(h2,ol) 


(h2,c2) 


s 


0.44 


0.31 0.36 


0.35 


0.21 


0.25 


0.27 


0.29 


Px 


0.02 


0.10 0.08 


0.10 


0.11 


0.10 


0.11 


0.11 


Py 


0.02 


0.10 0.09 


0.10 


0.11 


0.10 


0.11 


0.11 




0.14 


0.08 0.02 


0.03 


0.06 


0.07 


0.05 


0.06 


d2xy 


0.08 


0.04 0.07 


0.07 


0.05 


0.06 


0.04 


0.05 



(17) 



s 

Px 
Py 

d2xy 



Gi 



(hl,cl) (hl,c2) (h2,cl) (h2,c2) 

0.11 0.00 0.00 0.00 

0.00 0.00 0.00 0.00 

0.00 0.00 0.00 0.00 

0.00 0.00 0.04 0.00 

0.05 0.00 0.00 0.00 



G2 



(hl,cl) (hl,c2) (h2,cl) (h2,c2) 



0.15 0.01 0.11 0.02 

0.00 0.00 0.01 0.00(18) 

0.00 0.00 0.01 0.00 

0.00 0.00 0.00 0.00 

0.03 0.00 0.00 0.00 



Obviously, when it comes to the interband vertex (W) 
as well as all the other vertices, the "s-wave" chan- 
nel dominates, which retroactively justifies our earlier 
idealized analysis in terms of an attractive Hubbard 
model and a fictitious "superconductor." We find that 
- depending on the overall strength of Coulomb re- 
pulsion - the most likely ground state is a multiband 
VDW (SDW/CDW/ODW) which, due to the mismatch 
of the hole and electron bands and the underlying lat- 
tice effects, is generically in the "FFLO" region of the 
phase digram in Fig. [1] leaving portions of the origi- 
nal FS ungapped and metallic. As expected, this VDW 
(SDW/CDW/ODW) symmetry breaking at wavevector 
M is fueled by a large_susceptibility of nearly-nested hole 
and electron valleys [l8i| . 

If the Coulomb repulsion is just below what is needed 
to produce a metalhc VDW (Fig. [J), the mixed "Joseph- 
son" vertex G2 is strongly enhanced, as illustrated by 
(fTSj) and surrounding discussion. This is the regime 
where the interband superconductivity [1, [3, [i^] is pos- 
sible. Here, an important point needs to be made: there 
are two ways in which G2 can lead to high temperature 
superconductivity in Fe-pnictides: first, G2 itself can be 
the source of superconductivity. This is a naturally ap- 
pealing theoretical scenario, since it relies on the proxim- 
ity to a VDW (SDW) instability to enhance G2 and uses 
purely electronic interactions to generate superconduct- 
ing order. The difficulty in this case is that G2 has to 
overwhelm the intraband Coulomb repulsion and U'^ 
before superconductivity becomes possible, the condition 
being roug ly G2 > VTMj^. While G2 takes oflf as one 
approaches the VDW (SDW/CDW/ODW), there is a re- 
flection of this enhancement in the renormalized values of 
and as well. For the realistic model with four in- 
eqivalent bands and the interaction vertices displayed in 
Tables (|17l) and (fTS]) . this balancing act between J7's and 
G2 becomes very sensitive, particularly since the bare J7's 
start as generically larger and only two of G2's are not 
negligible in size while all four J7's are appreciable. Any 
effort to extend Eqs. fTS]) to four realistic bands and to 
all (intra, interband and mixed) vertices quickly descends 
into impenetrable numerics with the above sensitivity to 
the bare values in (fTTl - [T5|) , making it difficult to reach 
firm quantitative conclusions. A notable recent progress 
along these lines was made in Refs. [28j . 



All the vertices (fT7|) and ([T8|) are given in the original c, d 
electron basis ([TJS]) and are all positive (repulsive); they 
are easily converted into the basis used in pT|) by the p-h 
transformation (i.e., W —W, etc.). The numbers in 
the above table change if additional forms of interaction 
in real space are considered, for example those of Eq. ([4]), 
as discussed earlier. Again, provided one is outside the 
regime dominated by the Hund's coupling, such changes 
are minor. 



However, there is a reasonably straightforward way to 
illustrate the qualitative argument for the interband su- 
perconductivity mechanism near the VDW phase bound- 
ary. This argument follows straight from Eq. . Imag- 
ine that our isospin label is simply an ordinary spin. 
Therefore, we have spinful electrons with two (instead 
of four in real pnictides) orbital flavors, c and d [h and 
e). In this case, G2 is the interband pairing resonance in 
the spin-singlet channel. As argued earlier, we can safely 
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set Gi = and rewrite the remaining parts of p5)) as 



9u 
9w 
9w 

92 



-9l 

{g'wf 

{9'{v?. 



9I, 



= -2g2gc/ + 252 (5m/ + 9w) 



(19) 



where g's are functions of In(^) and g = dg/dhii^^^ and 
we again assume Uh = Ue = U . 

Imagine for the moment that there is no last term in 
the last hne of (fTO]) . i.e., W's and G2 are not directly cou- 
pled. As one moves to low energies w — s- 0, gw^s rapidly 
grow and one ultimately reaches the point where the sys- 
tem turns into a VDW. Meanwhile, U and G2 do nothing: 
this is easily seen by adding and subtracting the first and 
the last lines of (fT9|) gu ±g2 = — {gu ±52) ■ This is just 
the lowest order description of the real Cooper pairing 
instability in the s-wave spin-singlet channel, with the su- 
perconducting gap parameter having either the same sign 
for both bands c and d (conventional s-wave) or the oppo- 
site sign (s±-wave or an extended s-wave, s'), the corre- 
sponding coupling constants being gu+g2 and gu—g2, re- 
spectively. Both of these coupling constants are repulsive 
and thus both scale toward zero and into the Fermi liquid 
regime, leaving the VDW and W to determine the physics 
at low energies, unless G2 > \fThjTe) at the hare 

level. This is not impossible, but appears to be unlikely 
within the regime of interactions considered here, where 
G2 is typically quite a bit smaller than U . This tells us 
that the direct coupHng of G2 to VF's in the last line of 
(fTO)) must be crucial: the growth of W's as we approach 
the VDW eventually pulls G2 along with it, while U still 
continues being renormalized downward. This growth of 
G2 generated by its coupling to W's and the VDW could 
ultimately result in Gj > C/*(~ ^JU^U*), where Gj and 
U* are the renormalized coupling constants at some low 
energy scale ^ A, even though G2 < t/(~ \/UhUe) at 
the bare level. This implies that the coupling constant 
gu~g2 is attractive for oj < loo and translates into growth 
of s± (s') pairing correlations at yet lower energies. All 
of this is for nothing, however; Ws and the VDW insta- 
bilities are still far stronger. But, if the strong growth of 
W and the VDW instability are cut off by our fictitious 
"Zeeman splitting" in doped or pressurized FeAs (Fig. 
[1]) , then W stops growing at some energy scale ujz directly 
tied to the difference in FS size between h and e bands 
and the corresponding lack of perfect nesting. Then, if 
LJz < Wo, the subleading instability would take over and 
the ground state would be an s± (s') superconductor, ei- 
ther adjacent to the VDW boundary or coexisting with it 
in the pockets left ungapped by the VDW (Fig. [1]). The 
above argument is similar in spirit to the weak coupling 
mechanism for d-wave superconductivity once the single 
band repulsive Hubbard model is doped away from half- 
filling and the SDW ground state. There is an important 
difference, however: iiujz > wq there will be no supercon- 



ducting ground state since gu — 52 is still repulsive. This 
is a qualitative point and it underscores the fact that an 
s± (s') superconductor still has an overall s-wave symme- 
try and, unlike the nodal d-wave, must contend with the 
strength of the bare intraband repulsion. 

The second way is now obvious: G2 enhancement near 
the VDW instability can overcome the repulsion U if an 
attractive intraband interaction is at work as well. Such 
intraband attraction might come from phonons, for ex- 
ample. This attraction may or may not suffice to produce 
superconductivity by itself - the key point is that it re- 
duces the effective J7's PT)) allowing the enhanced G2 to 
cross over the hurdle. Note that in both of these cases, 
the purely electronic and the phonon-assisted one, the 
superconducting gap on the hole and the electron por- 
tions of the FS will have the opposite sign 0, [3, H^] , 
reflecting the fact that the interband pairing term G2 is 
repulsive. 



CONCLUSIONS 

In summary, we have proposed an idealized model of 
Fe-pnictides which includes an electron and a hole band, 
and takes advantage of their similar shape and size. The 
p-h transformation maps this model into a fictitious at- 
tractive Hubbard model in external "Zeeman" field. The 
ground states, fictitious superconductor and the "FFLO" 
state, correspond to insulating and metallic VDW in real 
materials. Next, by considering deviations from perfect 
nesting, two hole and two electron bands, and other re- 
alistic features of Fe-pnictides, we analyze the structure 
of interactions in the 8-1-8 orbital model 18| and identify 
the interband pair resonance mechanism that can gener- 
ate the real superconductivity in the region of the phase 
diagram of Fe-pnictides where the VDW order gives way 
to strong VDW fluctuations. 
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